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We present measurements of branching fractions of B° decays to multi-body final states containing 
protons, based on 232 million T(4S) — > BB decays collected with the BABAR detector at the 
SLAC PEP-II asymmetric-energy B factory. We measure the branching fractions B(B° — > D°pp) = 
(1.13 ± 0.06 ± 0.08) x 10~ 4 , B{B° -> D*°pp) = (1.01 ± 0.10 ± 0.09) x 10~ 4 , J3{B° -> D~ppTv+) = 
(3.38 ± 0.14 ± 0.29) x 10 -4 , and B(B° -> D*~pp-K + ) = (4.81 ± 0.22 ± 0.44) x 10" 4 where the 
first error is statistical and the second systematic. We present a search for the charmed pentaquark 
state, O c (3100) observed by HI and put limits on the branching fraction B(B° —* Q c pTY + ) x B(Q C — > 
D*~p) < 14 x 10~ 6 and B{B° -» e c p7r+) x B{Q C -» D~p) < 9 x 10" 6 . Upon investigation of the 
decay structure of the above four B° decay modes, we see an enhancement at low pp mass and 
deviations from phase-space in the Dp and Dp invariant mass spectra. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



The observations of the B — > D* pp7r + [lj and B — > 
D*-pn decays by CLEO 0, and the B° -> D°pp and 
B° — ► D*°pp decays by Belle suggest the dominance 
of multi-body final states in decays of B mesons into 
baryons Q compared to two-body decays. In this paper 
we present measurements of the branching fractions for 
the following four decay modes: B° — > D°pp, B° — > 
D*°pp, B a ->■ D-ppir+, and B° -> D*-ppir+. The study 
of the modes presented here can help clarify the dynamics 
of weak decays of B mesons involving baryons |5j ■ 

Since the branching fractions of multi-body decays 
are large it is natural to ask whether such final 
states are actually the products of intermediate two- 
body channels. If this is the case, then these initial 
two-body decays could involve proton-antiproton bound 
states (pp) 0, 01, or charmed pentaquarks 0, 0], or 
heavy charmed baryons. Motivated by these consider- 
ations, in particular the claim of a charmed pentaquark 
at 3.1 GeV/c 2 by the HI collaboration[ll|, the invariant 
mass spectrum of the proton-antiproton and the invariant 
mass spectra of the charmed meson and proton are inves- 
tigated. Throughout this paper, we shall use the terms 
"exotic" and "non-exotic" to refer to the "Dp" pair with 
total quark content cquud and cquud respectively (where 
q is u or d). Specifically, the "exotic" combinations refer 

to D ( *> p and D p while the "non-exotic" combina- 

, . — (*)o_ 

tions are D ( *> p and D p. 

The data used in this analysis were accumulated with 
the BABAR detector ^ at; the PEP-II asymmetric- 
energy e + e~ storage ring at SLAC. The data sample con- 
sists of an integrated luminosity of 212 ± 2 fb _1 collected 
at the T{4S) resonance corresponding to (232 ± 3) x 10 6 
BB pairs. The BABAR detector consists of a silicon 
vertex tracker (SVT) and a drift chamber (DCH) used 
for track and vertex reconstruction, an electromagnetic 
calorimeter (EMC) for detecting photons and electrons, 
a Cherenkov detector (DIRC) and an instrumented flux 
return (IFR) used for particle identification (PID). The 



efficiency of the selection criteria is determined with large 
samples of GEANT-based [l3| Monte Carlo (MC) simulated 
signal decays. 

We select D° decays to K + tt~ , K + n~ir°, and 
K + ir~ ir + Tr~ and D~ decays to K + ir~Tr~. We select 
D decays to D tt° and D*~ decays to I? 7r _ . The B 
candidates are reconstructed from D or D* candidates 
combined with a proton and an antiproton track and a 
pion track if appropriate. The D candidates are required 
to have a mass within ±3cr of the D meson mass, m5 
|l4|. The mass resolution, a(mu), ranges from 5.1 to 
13.0 MeV/c 2 for different D decay channels, the worst 
resolution corresponding to the mode with a it in the 
final state. The D* candidates are selected by requiring 
the mass difference AM = (mrj^—mo) to be within 3cr of 
the nominal value, AM, where a ~ 1.0 MeV/c 2 . Particle 
identification is required on the proton, antiproton, and 
pion from the B, and on the kaon from the D decay, us- 
ing combined information from the energy loss, dE/dx, 
in the SVT and the DCH and the Cherenkov angle in 
the DIRC. The proton identification efficiency is roughly 
90% with a mis-identification rate of less than 2%. To 
suppress backgrounds of all kinds, vertexing probability 
requirements are imposed on the D and B candidates. 
In order to reduce background from e + e~ — > qq events 
(where q is a u,d,s, or c quark), the cosine of the angle 
between the thrust axis of the B candidate and that of 
the rest of the event | cos (0bt)\ is required to be less than 
0.9 and the ratio of the second to the zeroth Fox- Wolfram 
moments is required to be less than 0.35. 

We select events in the region 5.2 GeV/c 2 < ities < 
5. 3 GeV/c 2 and \AE\ <0 .1 GeV, where m ES = 
\J (s/2 + p r • p b ) 2 /Ey — p 2 B (y/s is the total center-of- 
mass energy, p B is the B meson momentum and (Er, 
p T ) is the T(4S) 4-momentum, defined in the labora- 
tory frame), while AE — pr -Pb/^/s — {pr = {Er, 
Pr)> Pb = (Eb, Pb))- T ne selection is kept loose be- 
cause these two variables are used in a maximum like- 
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lihood fit to extract the signal and background yields 
simultaneously. If there is more than one B candi- 
date passing these criteria for an event, the candidate 
is chosen that minimizes x 2 = (mn — rhr)) /a(rri£)) 2 + 
(AM - AM) 2 /a{AM) 2 for the modes B° -> D*°pp 
and B° — > D*~ppn + , and the candidate that minimizes 
X 2 = {tud — fho) 2 / cr{mD) 2 for the modes B° — > D°pp 
and B° — ► D~pp-K + . 

The background for these modes comes from e + e~ — > 
gg events and from _B decays other than those under con- 
sideration. In both of these cases, the background comes 
from selecting random combinations of tracks and thus 
does not peak in either AE or to_bs. The one exception 
is in the case of B° — > D*°pp, where there is a possibility 
of events such as B° — > D ppir that peak at the B mass 
in rriES- However, since the ir° comes from the other B 
decay in the event, the AE distribution does not peak 
strongly in the signal region. 

We perform an unbinned extended maximum likeli- 
hood fit to extract the yields. The variables ities and 
AE are used as discriminating variables to separate sig- 
nal from background. The data sample is assumed to 
consist of two components: signal events and combina- 
torial background events due to random combinations 
of tracks from both qq and BB events. For the decay 
B° -> D*° pp, a peaking component is added to account 
for B° — > D ppir events. 



FIG. 1: Fit projections of thes for (clockwise from top-left) 
B° -> D°pp (D° -> K+tt-), B° -» D*°pp (if -> K + tt~), 
B° -> D*- P ptt + (S° -f K+ir~), and B° -» D~ppir + (D~ -> 
K + iy~tt~). The dashed line is the background contribution 
and the solid line is the background plus signal. 
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In addition, the signal is split into correctly recon- 
structed events (Class I) and mis-reconstructed events 
(Class II). The Class II events are signal events where 
one or more of the tracks from the signal B decay is 
lost and a track from the other B decay is included 
in the reconstruction. The fraction of Class II events 
is determined from MC and varies from nearly for 
B° -> D°pp -> K+ir-pp to almost 50% for B° -> 
D*~ppK + — * 7r~ K + Tr~7r ppTr + . 

In the maximum likelihood fit, each component is mod- 
eled by a probability density function (PDF) of the two 
variables ities and AE, 



V = V{m ES ,AE). 



(1) 



The likelihood for the N candidates in the event sample 
is given by: 



C = e- N '-H{N Big -[frV}+f ir r} I } 



-AWKkJ, (2) 



where N' is the sum of the fitted number of signal (iV s i g ) 
and background (iVbkg) events. The background PDF is 
given by 'Pbkg, Vi and Vn are the PDFs of Class I and 
II events in signal respectively, and // and /// are their 
corresponding fractions. 

The Class I signal events are parameterized with a dou- 
ble Gaussian for both m E s and AE. For Class II events, 
"niES is parameterized with the correlated function 
P H (m ES ,AE) = G(m ES )G 1 {AE) + P(m E s)G 2 (AE) 
where G represents a Gaussian and P a polynomial func- 
tion. All parameters for the signal PDFs are obtained 
from signal MC and fixed in the fit with the excep- 
tions of the means of the narrow components of the 
double- Gaussian distributions for both tues and AE for 
Class I events, which are allowed to vary. The com- 
binatorial background is parameterized with a thresh- 
old function 1 16| in vties and a second-order polynomial 
in AE, and all of the parameters are varied in the fit. 
The peaking background component coming from B de- 
cays in the B° — * D*°pp modes is modeled with a non- 
parametric 2-dimensional PDF in tues and AE and the 
yield is free in the fit. The tues distributions for the data 
and the fit, after selecting events with \AE\ < 20MeV, 

are shown in Figure □ for the D -> K+tt- and D~ -> 
K + tt~tt~ decays. 

For each event a signal weight is defined as follows: 



a 2 V 1 - 

sig sig 



cov(sig, bkg)P bkg 



(3) 



following the method described in Reference In 
Equation "Pg ig (^bkg) ^ s tne vaiu e of the signal (back- 
ground) PDF for event i; a s i g is the standard deviation of 
the signal yield; and cov(sig, bkg) denotes the covariance 
between JV s i g and iVbkg, as obtained from the fit. The 
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normalization of W* ig is such that their sum equals the 
total number of signal events, N sig . The sum of W* lg over 
a small area of phase space gives the correct distribution 
of signal in that area. 

The branching fraction is obtained as: 



where the sum is over all events i, N B g is the number of 
BB pairs in the sample, e$ is the efficiency for event i, 
which depends on its position in phase space, and B su h 
is the product of the branching fractions of the charmed 
meson decays 0,0]. We assume that the T(4S) decays 
with equal probability to B°B° and B + B~ . The statis- 
tical error on the branching fraction is obtained from the 
fractional error on the signal yield as calculated from the 
fit. 

The largest source of systematic error arises from the 
uncertainty in the charged track reconstruction efficiency 
determined from the MC. This systematic error ranges 
from 3.3% to 8.8% depending on the number of charged 
tracks in the decay mode. In addition there is a sys- 
tematic error due to the modeling of the PID efficiency 
for the protons and kaons of 4.5% for all modes and an 
additional error of 2% for the pion identification for the 
modes B° — > D~ppir + and B° — > D*~pp-K + . The un- 
certainty due to ignoring correlations between rriES and 
AE is estimated to be a few percent by performing fits 
to Monte Carlo samples that consist of fully simulated 
signal events embedded with parameterized background 
events. The uncertainties related to modeling of the sig- 
nal PDFs are calculated by allowing the AE and ities 
signal shape parameters for the £>° — > D~ppir + mode to 
vary in the fit and then varying the fixed parameters in 
the other modes by the differences observed between data 
and MC in this mode. This error ranges from 0.2% to 
2.8%. The fraction of Class II events is varied by 5% per 
7T°, or 5% for modes with no 7r°, to account for the uncer- 
tainty due to mis-reconstructed events and the difference 
observed is 1% to 5%. The uncertainty arising from bin- 
ning the efficiency in phase space gives a typical error of 
3%. Finally, the errors on the branching fractions of D 
and D* decays are included in the systematic uncertainty 
and range from 2.4% (for B° — > D°pp,D — » K + n~) to 
6.2% (for B° -> D* Q pp,lf -> K+n-it ). The total sys- 
tematic error ranges from 6.3% to 13.3%. 

The fitted signal yield and the measured branching 
fraction for each decay mode is given in Table [Q Av- 
eraging the branching fractions of the different D decays 
weighted by their errors and accounting for correlations, 



we obtain: 

B{B°^D°pp) = (1.13 ± 0.06 ± 0.08) x 10" 4 

B(B° -> D*°pp) = (1.01 ± 0.10 ± 0.09) x 10" 4 

B(B° -> D-ppn + ) = (3.38 ± 0.14 ± 0.29) x 10" 4 

B(B° -> D*-ppn + ) = (4.81 ± 0.22 ± 0.44) x 10" 4 

where the first error is statistical and the second system- 
atic. 

We investigate the decay dynamics by projecting the 
branching fractions obtained with Equation 01 onto the 
different invariant mass axes. This method requires that 
the variables used in the fit are uncorrelated to the vari- 
able being projected. The correlations between the in- 
variant masses and AE and rriES are observed to be 
small. Figure [21 shows the two dimensional projections 
(the Dalitz plots for the 3-body decays) for the four 
modes under study. Figure|3|shows 1-dimensional projec- 
tions and the comparison with phase space distributions 
for the pp, Dp (non-exotic minimal quark content of cud 
or edd) and Dp (exotic minimal quark content of cuuud 
or cduud) invariant masses. 

In comparison with phase space, an enhancement 
at low pp mass is seen in all decay channels. Such 
an enhancement has been observed in other situations 

EM E3, H3, 

indeed, it is also observed in the back- 
ground pp distributions in this analysis. In the left plot 
of Figure the pp distributions for all four modes have 
been overlaid removing the events with M(D p) less 
than 3.1 GeV/c 2 and normalizing to the total area. In 
addition, each event entering Figure^has been weighted 
by a phase-space factor and thus the distribution is pro- 
portional to the square of the matrix element. The dis- 
tributions of the four modes show the same behavior. We 
have also compared our phase-space corrected pp distri- 
butions (averaged over the four modes) to those measured 
in e+e" -> ppj Hj and B+ -> ppK+ HH by BABAR, 
shown on the right of Figure 21 arid again there appears 
to be good agreement. 



TABLE I: The branching fractions (in units of 10~ 4 ) for the 
B° decays considered here. The first error is statistical and 
the second systematic. 



B" decay 


D decay 




23(l(n 4 ) 


B° -» D°pp, 




214±16 
514±38 
320±26 


1.09±0.08±0.08 
1.15±0.08±0.10 
1.24±0.10±0.11 


B° D*°pp 
D*° - D°n° 


K+tt- 
K+n-n° 


57±9 
104±19 
46±12 


1.21±0.17±0.11 
1.08±0.14±0.14 
0.75±0.18±0.09 


B" -» D-ppn + 


K + 1T~%~ 


1166±47 


3.38±0.14±0.29 


B° -» D*-pp7r+, 
D*~ -> D a %- 


K+TY~TT° 


241±18 
522±32 
311±24 


4.84±0.40±0.44 
4.71±0.30±0.50 
5.05±0.42±0.59 
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FIG. 2: The distributions of the branching fractions (in units of 10" 6 / GeV/c 2 ) for B° -» D°pp (top row, left), B° -> D*°pp (top 
row, rig ht), B° -> IT "pp7r (middle row), and B° -> D*~ ~ppn (bottom row) projected over two invariant mass dimensions. 
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Explanations that have been proposed to account for 
the enhancement observed at the pp threshold include 
a gluonic resonance 1231 a nd short-range correlations be- 
tween the p and the p |24| . The BES collaboration has re- 
cently claimed evidence for a resonance decaying to ttttt]' 
with a mass of 1834 MeV/c 2 and a width of 69 MeV/c 2 
[1^1 . This resonance should also decay to pp and the mass 
and width measured by BES in imr]' is in agreement with 
the enhancement seen by BES in the pp distribution in 
3/tp — > jpp decays |22j| assuming a Breit-Wigner with 
corrections for final state interactions [2H l2?j | . 

With respect to the Dp invariant mass spectra, other 
than an excess at low mass in the B — > D°pp mode, 
the plots in the middle row of Figure [3] are in qualita- 
tive agreement with the phase space histograms. The 
low mass excess in B° — ► D°pp is also easily seen in 
the Dalitz plot in Figure and appears again to be a 
threshold enhancement as in the pp case. While it would 
be expected that the same effect would be seen in the 
B° — > D*°pp mode, the statistics are much lower and the 
mass threshold is higher. 

The Dp distributions, in the bottom row of Figure 

— (*)o 

we observe a clear tendency to peak toward high D p 
mass in comparison with phase space for the three-body 



modes. This is also reflected in the apparent asymmetry 
in the Dalitz plots. The four body modes are in qualita- 
tive agreement with phase space distributions in the Dp 
projections. 

The HI Collaboration has claimed evidence for a 
charmed pentaquark state decaying to D*~p at 3.1 
GeV/c 2 whose width is less than their experimental reso- 
lution of 7.1 MeV/c 2 . By fitting the D~p invariant mass 
spectrum in the decay B° — > D~ppn + to a Breit-Wigner 
plus linear background, we obtain an upper limit on the 
branching fraction: 

B{B° -> e c jJ7r+) x B(Q C -> D~p) < 9 x 1(T 6 , (5) 

while for the D*~p spectrum in B° — > D*~ppTT + we ob- 
tain: 

B(B° -» e c p7r+) x B{Q C -> D*-p) < 14 x 1(T 6 (6) 

at 90% C.L. For this limit we have assumed the resonance 
width for the <d c to be 25 MeV/c 2 , which corresponds to 
the upper limit on the width given by HI. If we assume 
a smaller width, the limits decrease. 

In conclusion, we have measured the branching frac- 
tions of B° -> D°pp, B° -> D*°pp, B° -> D-ppTi + , and 
B° — > D*~ppir + . The results obtained for the modes 
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B a _^ D*-ppn+, B° -> D*°pp, and B° -> D>jj agree 
with the previous measurements and have smaller uncer- 
tainties while the decay B° — > D~ppir + has been mea- 
sured for the first time. We do not observe any evidence 
for the charmed pentaquark observed by HI at M{D*~p) 
of 3.1 GcV/c 2 . In comparison with phase space we ob- 
serve a low-mass pp enhancement similar to other obser- 
vations in pp production. We also observe a deviation 
from phase-space structure in the Dp and Dp invariant 
mass distributions for the three-body modes. 
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FIG. 3: The branching fraction (B, in units of 10 6 / GeV/c 2 ) distributions versus pp (top), non-exotic (i.e Dp) (middle), exotic 



(i.e Dp) (bottom) invariant mass for (from left) B" -> D"pp, B" -> D*"pp, B" -» D"pp7r + , and B° 
decay modes combined. The solid lines are the distributions expected from a purely phase-space decay. 

^20 = 
-100- 



D* 



ppn + with all D 




2 2.2 2.4 2.6 2.8 3 3.2 3.4 
m(pp) (GeV/c 2 ) 



2 2.2 2.4 2.6 2.8 3 3.2 3.4 
m(pp) (GeV/c 2 ) 



2 2.2 2.4 2.6 2.8 3 3.2 3.4 

m(pp) (GeV/c 2 ) 



2.2 2.4 2.6 2.8 3 3.2 3.4 
m(pp) (GeV/c 2 ) 



n 



-1 













14^ 
12^ 
10r 

8^ 
6r 
4r 



3 0- 

« -2L 




3 3.2 3.4 3.6 3.8 4 4.2 
m(D°p) (GeV/c 2 ) 



3 3.2 3.4 3.6 3.8 4 4.2 
m(D*°p) (GeV/c 2 ) 




2.8 3 3.2 3.4 3.6 3.8 4 4.2 

m(D p) (GeV/c 2 ) 



3 3.2 3.4 3.6 3.8 4 4.2 
m(D* p) (GeV/c 2 ) 




18? 
16^ 
14? 
12 L 
10? 

8r 

6? 
4? 



ffl 0- 




3.2 3.4 3.6 3.i 

m(D°p) (GeV/c 2 ) 



3 3.2 3.4 3.6 3.8 4 4.2 
m(D*°p) (GeV/c 2 ) 



2.8 3 3.2 3.4 3.6 3.8 4 4.2 
m(D"p) (GeV/c 2 ) 



2.8 3 3.2 3.4 3.6 3.8 4 4.2 

m(D*"p) (GeV/c 2 ) 



FIG. 4: Left: The phase space-corrected pp invariant mass distributions for all four decay modes: B° — > D°pp (triangles), 
B° °PP (open circles), B° — » D ppn + (squares), and B° — > D* pp7r + (closed circles). Right: The pp distributions from 

the present analysis averaged over the four decay modes (closed circles) compared to the distributions obtained in e + e~ — > ppy 
(open squares) and B + — > ppK + (open circles). 
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